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We continue our study of the statistical properties of particles in equilibrium obey- 
ing Smoluchowski dynamics. We show that the system is governed by a kinetic 
equation of the memory function form and that the memory function is given by 
one of the self-energies available via perturbation theory as introduced in previous 
work. We determine the memory function explicitly to second-order in an expansion 
in a pseudo-potential. The method we use allows for a straightforward computation 
of corrections via a formal expansion and we therefore view it as an improvement 
over the conventional mode-coupling theory (MCT) formalism where it is not clear 
how to make systematic corrections. In addition, the formalism we have introduced 
is flexible enough to allow for a wide array of different approximation schemes, in- 
cluding density expansions. The convergence criteria for our formal series are not 
worked out here, but the second order equation that we derive is promising in the 
sense that it leads to analytic and numerical results consistent with expectations 
from computer simulations of the hard sphere system in addition to replicating the 
desired features from conventional MCT (e.g., a two-step decay). These particular 
solutions will be discussed in forthcoming work. 

I. INTRODUCTION 

A powerful approach for studying the dynamics of systems of classical particles was pre- 
sented in Ref. 1 (referred to here as FTSPD). In Ref. (referred to here as SDENE) 
this method was developed to study fluctuations in equilibrium for systems obeying Smolu- 
chowski dynamics (SD)^. A parallel development is carried out in Ref. for systems obeying 
Newtonian dynamics. The treatment of SD is extended here to include a full self-consistent^ 



2 



treatment of the separation of statics and dynamics and a derivation of a kinetic equation 
- valid to second order in perturbation theory - governing the density fluctuations over the 
entire time range. 

In SDENE, the fluctuation kinetics are described in terms of a self-energy structure. 
These self-energies can be conveniently obtained in a perturbation theory expansion in a 
pseudo-potential. Furthermore, the self-energies are divided into single-particle and collec- 
tive contributions. In SDENE, we focused on the collective contribution which governs the 
long-time slow kinetics in the problem, and showed that one can find a simple self-consistent 
relationship between the static structure factor and the zero-frequency component of the 
collective part of the self-energies. 

In this paper, we fully analyze the single-particle contributions to the self-energies. We 
show the following: 

1. The single-particle contribution to the self-energy can be associated with the equation 
of state governing the system. 

2. While the collective degrees of freedom dominate the long-time dynamics in this sys- 
tem, the single-particle degrees of freedom govern the early-time kinetics and the 
approach to the slow-regime. 

3. We derive here the kinetic equation of the memory-function type valid at second-order 
in perturbation theory and including both single-particle and collective contributions. 

It is our intention to solve this kinetic equation numerically in future work^, and to 
show that the analysis of the collective contribution in SDENE can be extended analytically 
to obtain a two-step kinetic process similar to that obtained from mode coupling theory 
(MCT}i. 

MCT represents the current de facto theoretic description of dense fluids and the transi- 
tion from fluid to glassy state^"— . However, MCT is limited by its ad hoc construction and 
lacks a mechanism to institute systematic corrections. We derive here, on the other hand, 
a form of the kinetic equation of the memory-function type used in MCT, and our theory 
provides the crucial advantage of well-defined, perturbative corrections. Thus, one can use 
our methods to study vertex corrections, three- and higher-mode corrections to the standard 
MCT two-mode form, and high-frequency effects. 
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We begin with a brief review of the previous work to estabhsh context and notation. 
We next develop the equation of state described above and then complete the second-order 
vertex function. In the final section, we derive the kinetic equation and discuss the memory 
function at its heart. 

II. REVIEW OF PREVIOUS WORK 

The key components of the field theory approach are the two-point matrix cumulant 
functions, Gij{q,uj), and the two-point irreducible vertex functions, rij{q,uj), with i and j 
running over the fields p and B where p is the particle density and S is a response field. 
The kinetic equation of interest results from an analysis of Dyson's equation which takes 
the form 

^ikGkj = Sij. (1) 

k 

In SDENE it was shown that the two-point irreducible vertex can be separated into two 
contributions, 

Tij = jij + Kij, (2) 
where is the single-particle contribution and K^j is the collective contribution. The 

(2) (2) 

second-order contribution, K^j , was derived and partially analyzed. It was shown that K>j ' 
itself satisfies a fluctuation-dissipation relation while remaining a quadratic functional of the 
two-point matrix correlation function, Gij. The single-particle contribution to the two-point 
vertex, 7jj, is defined as the inverse of the quantity 

G,^ = Tr 0,0,6^-^+^^ (3) 

such that 

likGkj = Sij (4) 

k 

where the fields $ = (p, B) are one-particle additive, 

N 



a=l 
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H = {Hff,, Hb) is a conjugate external coupling field, and the term /SW contains the pseudo- 
potential interaction. 

The form of the interaction is defined carefully in FTSPD, but to second order in 
the pseudo-potential is given by 

AW = AW^^^ + Aiy(2) + (6) 

where 

Al^(^) = (7) 

u 

Aiy(2) ^ 1 (8) 

u,v 

and 

Fi = ^crij(l)j, (9) 
j 

and where the interaction matrix aij is defined by 

(^ijiQ) = V{q){Sip5jB + SisSjp) (10) 

where V{q) is the Fourier transform of the potential. (We have not found it confusing to 
use the same symbol for both the coordinate- and wavenumber-space representations.) 

The trace, Tr, is defined carefully in FTSPD and SDENE, however we do not need the 
details in defining Eq.dH]) since all the noninteracting cumulants among the fields p and B 
are available for the non-interacting case in FTSPD. This suffices to determine Qij,,,k- We 
gave the solution for the noninteracting 7!°'' in both FTSPD and SDENE and provide a 
summary of the results in Appendix |X1 

In general, the procedure is to determine 7jj and Kij, then solve Eq.([T]) to obtain the 
two-point correlation functions Gpp, GpB-, and Gsp- In the special - but very important - 
case where the system is in thermal equilibrium, one finds that there is a simple fluctuation- 
dissipation relation (FDR) between the two-point quantities, 

Gpsiq, w) - Gspiq, w) = if3uGpp{q, u). (11) 

In SDENE, we showed that one has in this case a simple kinetic equation satisfied by the 
density- density correlation function 

^Gpp{q,t) = -Dq'pS-\q)Gpp{q,t) - Dq' j^' ds/3'pEBB{q,t - s)^Gpp{s) (12) 
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where T.BBi.lit) is the kinetic contribution to occurring in Eq.(IT]), S{q) is the static 
structure factor, p is the average density, and D = ksTD is the product of the temperature 
and diffusion coefficient. In the same work, we determined the collective contribution to 
Vij to second-order in an expansion in a pseudo-potential and established that F^^ and 
themselves satisfy a FDR and are quadratic functionals of the exact density-density 
correlation function Gpp{q,t). Here we want to determine the "single-particle" contribution 
to Tbb- 



III. EQUATION OF STATE 

In our approach here, we generate approximations for both the static and dynamic prop- 
erties. In SDENE, we showed how approximations for the static structure factor entered the 
analysis and we used the collective part of the self-energy, Kij, to make contact with the 
equilibrium statics via the static structure factor. Here, we want to show how the equation 
of state enters the development. 

A. Equation of state 

In FTSPD, we established the fundamental identity for the one point quantity 

Gi = Tr 0^6^-^+^^ (13) 

where i labels space, time, and fields p or B. This is the equation of state. For zero external 
field, H = 0, and keeping terms to second-order in the pseudo-potential as given by Eq.Q, 
we find the one-point quantity 



Gi = Tr (pi 



Tr ^i + TT4>iJ2 FuGu + Tr <Pi\ Yl FuF,{Guv + GuG,) + 



2 

2 



u,v 



Tr 0i + ^Tr (piO'uk<PkO've(peiGuv + GuGv) + ... 



2 

u,k u,v,k,£ 



{Gu, + G^G,) + ... (14) 



2 

u,k u,v,k,£ 
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Using wavenumber and frequency labels such that 1 = {ki,ui), we have the zeroth-order 
contribution (valid for uniform systems) 

Gf\q,,u,) = 5,,5il)po (15) 

where we introduce the notation 6{1) = {2iTY6{qi)27[6{uJi) and where po is the density in 
the absence of interactions in the grand canonical ensemble. 
Next, we have the first-order contribution 

G['\l) = Glf{12)a,,{23)G,{3) (16) 

where we now move to a convention where summation over repeated indices and integration 
over repeated, barred variables are implied. At all orders, the two point cumulant has the 
form 

G.,(12) = G,:,(l)5(l + 2) (17) 

(with 6{1 + 2) = (27r)°'5(gi + q2)2TT6{ui + 002)) due to translational invariance, and likewise 
the full one-point is 

G,(l) = 6,,6{l)p (18) 

where p is the average density. 

It is easy to show that the first-order contribution then yields 

^^^(l) = GS(12)aBp(23)G,(3) 
= gS(10)F(0)G,(0) 
= {-PGf\l))V{0)-p 

= -f3pV{0)poS{mp (19) 

where, in the next to last line, we used the identity G^^^{10) = —l3Gp{l) discussed in 
Appendix [Bl 

Turning to the second-order contribution, we have 

G<p)(i) = G'f'^)(l) + G'f'^)(l) (20) 



where 



Gf''\l) = ^G;°i(123)a,.(24)a,.(35)G.(4)G,(5) (21) 
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and 




-Gi;^(123)or,„(24)afc,(35)G'™(45). 



(22) 



We can express the three-point cumulants in terms of the three-point irreducible vertex, 



(The non-interacting three-point vertex functions 'Jljl are summarized in Appendix [XI) This 
gives, 



^^'■'^(l) = -^GS)(16)Gg)(27)Gl°i(38)7SU678)a,.(24^ 
= -^gS)(16)7SU678)G«(7)G«(8) 
= -^G£)(16)7SU678)(-/3V(0)p)Gf (7)(-/3V^(0)p)Gf(8^ 
= -^GS)(16)(/3\/(0)/.)2(-7(J(67)Gf(7)) 
= ^(/3V(0)p)2Gf (1) 

= lmiO)pfpo5.,5il) (24) 



where (in the fourth hne) we have used the identity (estabhshed in Appendix [B]) 



Gi;;i(123) = -G,^^^(14)Gi";(25)Gl"j(36)7SU456). 



(23) 



7SUl23)Gf(3) = -7i?(12). 



(25) 



Next, we have the more comphcated second-order contribution 



Gf''\l) = -Ggi(123)a,„(24)a,,(35)G„,(45) 
= -^g£^(12)7SI(234)G,.(34) 



(26) 



where we again have an effective propagator 



G,,(12) = Gi?(13)a,p(34)Gp„(45)a,,(56)G'l°)(62) = G,,(l)5(l + 2). 



(27) 



Looking at the i = B component, we have 



Gf'\l) = -2G21(12)7Sl(234)G,.(34) 
= -^G'gJ(12)7Sl(234)G,.(34). 



(28) 
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Enforcing the (5-functions, we find 



1 



d3 



7^°U0, 3, -3)G,5(3) + 7X(0> 3, -3)Gbp(3) 



(0) 



(29) 



Consulting the forms given in Appendix |A] for the three-point vertex functions, we see 
that each contributes at most a term hnear in frequency. When we perform the frequency 
integral over the response function weighted by this factor, we find that the integrals vanish 
since one can close in the half-plane where the response function is analytic. (This is the 
upper half-plane for GpB and the lower half-plane for Gsp-) Thus, 



G 



(2,2) 
B ' 



Returning to the p-component of Eq. 



1) = 0. 
, we have 



G?''\l) = -2G(°i(12)7SU234)G,.(34). 



(30) 



(31) 



The integrals over the response components of G vanish for the same reasons given above 
and one is left with the contribution 

1 



2<i(10) J rf37),^;^(0,3,-3)G,,(3) 
-^(-/3G(°)(1)) I rf37it(0,3,-3)G,,(3). 



„{o) 



Then, using 



we are left with 



7l?;,(0,3,-3) = l//3p; 



0' 



(32) 



(33) 



Gr\l] 



5(1) 



S{1)— I d3G,,{3) 



d'^k du ~ 



5(l)po/3 



2po J i^ny 271 
d'^k 



Gpp{k,uj) 



(2vr)^ 



V\k)Sik) 



where we use the result 



^Gpp{k,co)=plP'V'{k)S(k) 



(34) 



(35) 
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derived in Appendix O 

Collecting all the contributions to the equation of state, - Eqs. f|T^ . f lT9|) . (^^, and 
- we have 



2''^'^ ' 27 (27r)' 

= Po (1 - V^(O) + ^V^^(O) + 1 1 (36) 

where V = p(3V and S = S/p. This agrees with a strictly static formulation of the problen>i^, 
and can be rewritten in the form 

1 r d'^k ~. 



P 

— = exp 
Po 



(37) 



Expanding in V leads back to Eq.(l36ll. 

We elaborate on the equation of state and connect it to a more conventional form in 
Appendix [El 



B. Treatment of Qij 

We next need to treat the propagator 

g,, = Tr (/,,<^,e^-<^+^^[^] (38) 

in order to determine jij using Eq.dl]). We shall see that Qij is, roughly speaking, a single- 
particle quantity. The expansion in powers of V follows closely the expansion in treating 
the equations of state. After expanding AW in powers of V, we have 

g,, = + Gfla^tGt + \G%,a^^a,,{Gu. + G^G^) + 0{V''). (39) 

The first term is just the noninteracting matrix propagator determined in FTSPD as 

GSf(12) = GSf(l)5(l + 2) (40) 

with 

G»^(gS(1))^^ (41) 



and 



GS'(l) ^ (42) 



The first- and second-order terms are addressed in turn. 
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1. First- order 

For the first-order contribution, 

<3\f = G\%a,eGe, (43) 

we again replace the zeroth-order three-point cumulant with the three-point vertex using 
Eq.(l23]) and find 

= G^{-^^)G^ (44) 
where the self-energy is given in Fourier-space by, 

^11^12) = 7il(123)Gi°l(34)a,,(45)G,(5) 
= 7il(123)GS(30)G,(0)\/(0) 
= 7il(123)(-/3GW(3)p)\/(0) 
= -/3l^(0)p(-7i°Hl2)) 

= (3pV{0h^J{12). (45) 
Putting this back into Eq. lHil) and remembering that G^^juJc'^^j = Of^, we have 

g^^\q,Uj) = -Gf^{q,uj)pVm (46) 
Combined with the zeroth-order result, we have 

Qf^ ^) + ^i;^ (g, ^) = ^) [1 - ^(0)/5p] (47) 

and we see that the first-order contribution is a static contribution to the equation of state. 

2. Second- order 



Working at second-order we have the two-pieces contributing to 



6;f = 6w> + ar' (48) 



where 



^ir^ = \G%auucy,.G^G, (49) 
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and 



^(2,2) ^ ryiO) ^'K^^ 



Let us take the disconnected piece first. Using the representation = p5(l)5up, we 

have 

^ir^(12) = ^G,%(1200)(\/(0)p)^ 

= \ii''Gf^{\mmf (51) 

where we use the identity from Appendix [B] given by 

G!°L(1200)=/32Gg)(12). (52) 

Turning to the substantial contribution, we have 

^(2,2) 1^(0) it.o\ 

For the first time we encounter the noninteracting four-point cumulant, 

^(0) _ ^(0)^(0)^(0)^(0) (0) /r.N 
'^ijki ~ '^ix ^jy ^kz ^Iw Ixyzwy K"^^) 

where we have introduced the four-point vertex 'ylhjzw- It is important to recognize that 
the four-point vertex has a one-particle reducible contribution and a one-particle irreducible 
contribution, 

(0) ^ (o,R) ^ (0,1) _ /55^ 

ixyzw ixyzw ' Ixyzw \ / 

The reducible contribution can be written quite generally as 

lxy^2 = IxykGkafl + ifzi^kaSjl + Ixlk^kelill (56) 

where, to the order we are considering here, we may take Qij — )■ G^-^'' such that 

^{0,R) _ ^(0) ^(0) (0) (0) ^(0) (0) (0) ^(0) (0) 

Ixyzw — Ixyk^kl hzw ' Ixzk^kl hyw ' Ixwk^kl Uyz' W J 



The non-interacting four-point irreducible vertex functions are determined in Ref. 
summarized in Appendix 1X1 
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and 
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Using Eqs. (l55i) and (157|) in Eq. (lMI) . we obtain the second-order contributions to Qij, 

^(2,2) ^ ^(2,2,^,1) ^ ^(2,2,K,2) ^ ^(2,2,K,3) ^ ^(2,2,/) ^^g^ 

where we have four pieces. The first is a disconnected contribution 

^(2,2,K,l) _ 1(^.(0)^(0) (0)^(0) (0) ^ /c.q\ 

where G is defined by Eq. fl27p . The next two terms are one-loop contributions given by 

^(2,2,i?,2) _ 1^(0) (0)^ ^{0) (0) ^{Q) ,r,n.s 
^ij ~ 2 «^ Ixze'^zw'^es Iwsy^yj ) 

and Q^j'"^'^'^^ = Qij''^'^''^\ Finally, we have a Hartree-like contribution 

^(2,2,/) _ 1^(0)^(0) (0,/)^ ((^-.s 
^ij ~ 2 if Ixyzw'-^zw- K'^^J 

Consider first the disconnected term 

= ^G£)(13)GfJ(24)7Sk345)Gi?(56)7(l(678)G^^ (62) 
Since the frequency integral over the response components of Gzw vanish we have 

gif''''''\l2) = ^G£)(13)GgH24)7S),(345)G(?(56)7S 

= ^G!?(13)Gg)(24)7S),(345)Gfi(56)7l,t( 

= ^G;?(13)Gg)(24)7SU345)<i(56)7it(678)G^^ 

Let us pause and look at the combination G'^'^(56)7^j^(678)Gpp(78). Using the implicit 
delta functions, we can write this as 

<i(56)7Sp(678)G,,(78) = <i(5,0) J rf77i°jp(0, 7, -7)G,,(7) (63) 

Because 7spp(0, 7, — 7) = the only remaining integration is over (5pp(7). Inserting 

this, along with G'^'^(5,0) = —l3Gf\5), we have 

<i(56)7il(678)Gpp(78) = (-/3GW(5)) J d7G,,{7)/Ppl (64) 

which gives 

'^'^^(12) = ^G;?(13)G5y(24)7i5p(345)(-GW(5)) J d7G,,{7)/pl. (65) 
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If we finally use 7iyp(345)Gp°^(5) = — 7i°-*(34), then we are left with 

= loifiU) J d7G,,{7)/pl. (66) 

This term goes into the statics and the determination of the equation of state to second- 
order. Combining this with the zeroth- and first-order contributions, we have the very simple 
result 



MO) I I (-(2,1) , M2,2,R,1) ^ ^(0) 



2" ' ^ " 2 



= G^^. (67) 
Po 

The remaining three terms in Eq. flSHl) are of the form of self-energy terms which make 
contributions to 'flf. 



3. Summary of results at second-order 

The results for the single-particle propagator to second-order in perturbation theory can 
be written in the form 

a« = fG™-G.'M?Gg'. (68) 

Po 

where 

7?^ =7r"^+7r^ (69) 

and where 

(2,loop) _ (0)^ (0) ,r,n.-s 



and 



Ti, = -ifliMGu. (71) 



To this order in perturbation theory we can rewrite Eq.( l68|) in the form 

6. = - Gl^g'a,, (72) 
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where gI'^^ is Cf^^ with po replaced by p and 7^^'' is 7^^^ with po replaced by p. Comparing 
with Eq.dl]), we can identify 

7,= 7;°' +7.?'. (73) 
There is no explicit first-order term for the single-particle two-point vertex. 

C. The full two-point vertex function 

The two-point vertex is then the sum of the single-particle contribution 'jij and the 
collective contribution Kij, 

rij=iij + Kij. (74) 

Let us review the collective contributions and discuss the full vertex function. 
The first order collective contribution is simply 

Kff = (75) 

while the details of the second order collective contribution are worked out carefully in 
SDENE and result in 

■^in '^^iuv^uw^vz^ wzj (76) 



2 



where 



In SDENE, we kept terms of the lowest order in the vertices and effective cumulants to 
develop the nontrivial approximation 

with the symmetrized propagator 



^ij — iji^fk ^kiG ij + GikCkiGf^). (79) 



Bringing all these pieces together, we have for the second-order two-point vertex, 

r.,(i) = -if {I) + i^Hl) + 7{?^(l) + Kf{^- (80) 
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Notice that 7^^'''"^^ and K^j^ share the same one-loop structure, but with different prop- 
agators. All the propagators - Gij, G\'^\ Gij, and Gij - satisfy the fluctuation-dissipation 
relation. 

It was shown in SDENE that loop contributions like these can be written as 



and 



where 



Kg) ^ d^^^\j,j] (82) 

X J ^^G{q,,oo,)G{q,,u,)J,„ (83) 

x/^^G(g3,u;3)G(°)(?4,a;4)J,„ (84) 

J _ 1 ^1 (1 + ^(^3^13 + ^4^14))^ 



and 



UJ^-\- Ul4^ — Uli — IT] 

Let us look at the static hmit where the structure factor is related to the potential in 

perturbation theory by 

The terms up to first order are easy to simplify in the a; — > limit and we have 

"rPpM^-j-p (88) 

and 

K^liq.O) = -V{q). (89) 
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Looking next at the second order terms, we have for the collective loop term 



d'^'k?. duj?. ^ 



G'pp(g3,W3) 



dojA ^ 



2tx 



2/3p4 J (2Ti)d 2ti 

2^-4 J 0^,{-pPV{h)S{h)){-pPV{q - h)S{q - A:3)) 



7<1 



-0^V{k)~S{k)V{q-k)~S{q-k) 



(90) 



where we have used the zero time results derived in SDENE. 

Doing the same with the single-particle loop contribution, we find 



(2,/oop) 
IBp 



1 



(g,o) = - 

P 



d'^k^ 

d'^k^ 



^Gpp(A;3,W3) 



^^G^^iq-k„u,) 



d% 



/3p2 J (27r)' 



2n - 
V\q)~S{q) 



Gpp{k^,uj^) 



(91) 



We now look at the Hartree-like term. Using the results summarized in Appendix [XJ we 
have 

2 



^Bppp'^PP + ^BpBp'-'Bp + 7b 



(0,/) 
ppB^pB 



d'^k^ dus 



2 + 



iuja2 ^1^303 ^1^304 



2/3p3 J (27r)'^ 27r [v «^3 
— — G'Bp(/i;3,u;3) - (5^5(^3, W3) 

/3 Kg /3 ^3 



+ 



«3 y 



Gpp(fc3,a;3) 



where (under our constraints) 



(92) 

(93) 
(94) 
(95) 

Noting that integrals over uj^Gpp{uj^) vanish due to odd symmetry and integrals over the 
response functions vanish when the contour is closed in the appropriate half plane, we are 
left with 



^2 


= K^- 


[-3,1-^1- 


h3 -1 + -^1-3,1-^1-3,-1 


"3 


= Ki^ 


-3,1^1^ 


-3,3; 


0:4 


= Ki- 


h3,l-^l- 


h3,-3- 



iBp (9.^) 



d^k. ~ 



+ 



/3p2 J (2'n)d 
d'^k- 



V\k,)S{k,) 



lUJ 





2/3p2 J (2ttY\ki) 



V\k3)S{k,) 



(96) 



17 



We can take a closer look at the second term. Writing out the full form of q;2, we have 
-iuj 



d% f{,^ + ,- hf ^ (f-^- M^ ) v^^k,)Sik,) (97) 



[ d% f 
J 



If we perform a change of variables, 



p'^ ^ [q- kf = q^ + e-2q-k (98) 
^;2 = (g + kf = q'^ + k'^ + 2q- k, (99) 



then the measures of the two terms change to 

d^ks = dk^d(t)dukl = -'^^^dksdp = '^^^dksdv (100) 

where u — cos{9) and du — —sm{9)d9 is the standard angular integration substitution. 
Using this set of variables, it is easy to see that the term vanishes, 

d'^ks f Ci2\-rr2n \an \ f ^^3-fy2n \ an \ f f dv k-^v {v"^ + q^ - k'^Y 



^y\h)~s{k,)^ j'^nmk)[- j 



dp k^p (p^ + q'^ — k"^^^ 
2% q Ap^q 



+ / T-^ '' ' ) = 0. (101) 



This leaves us with an o'-independent Hartree term given by 

'<'BfM = ~j^,nk^)S^k,^ (102) 
which cancels the loop term: 

7g(?,0)=0. (103) 
Collecting the results, we have finally 



or, inverting, 

Jd 



~S-\q) = 1 + V{q) - ^ y ^,nk)S{k)V{<l " k)S{q - k). (105) 
This is the same quantity evaluated in SDENE to determine the effective potential. 
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IV. THE KINETIC EQUATION 

A. The kinetic equation 

Having determined Tij to second-order, we could proceed to solve the Dyson's equation 
for Gai3- However, there is a more economical route sketched out in SDENE that takes 
advantage of the FDR and which is nonperturbative. Using the FDR, we can go from 
Dyson's equations to a single equation for Gpp{k,t). We fill in the details of the analysis 
given in SDENE. 

To derive the kinetic equation, we begin with the Bp component of Dyson's equation: 

^BxGxp = ^bbGbp + ^ BpGpp = Sbp = 0. (106) 
In q-, t-space, this is explicitly 

j ds Tspiq, t - s)Gpp{q, s - t') + J ds TBsiq, t - s)GBpis - t') = 0. (107) 
Let us split the two-point vertex into two contributions as 

r,,{q,t)=^if{q,t)-E,,{q,t), (108) 
where we define 7^^^ (g, t) to be all terms local in time such that 

7ll\q,t- s)G,M^ s-t')= ^ll\q,t)Gk,iq,t- t'). (109) 
From this, it follows that 

The second group, Sjj(g, t), is the dynamic memory function^ which retains its convolution 
form and is made up of the remaining contributions, 

s.,(g,t) = -7r'''^-^^ (111) 

We may now write Eq. (11071) as 

^^Blit)Gpp{t,t')+l^Bl,{t)GBpit,t') = ^{t,t') (112) 

where 

^{t, t') = I dsJ:Bp{t - s)Gpp{s -t')+ I dsJ:BB{t - s)Gbp{s - t') (113) 

J —oo J —oo 
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using the fact that T^spit — s) ~ 9(t — s) and Gbp{s — f) ~ 9(t' — s). We then use the 
fluctuation-dissipation relations 



and 



j:Bp{t-s) = e{t-s)/3-^j:BBit-s) 



GBpis - t') = e{t' - s)l3^^Gpp{s - 1' 



(114) 



(115) 



to obtain 



ds 



d 

(3 — T.BB{t - s) 



Gpp{s -t')- J ds^BBit - s)P-^^Gpp{s - 0(116) 
If we integrate the first integral by parts, we have 

/t Q 
dsT.BB{t-s)—Gpp{s-t') 

dsT.BB{t-s)—Gpp{s-t') 

= -PT^BB^mppit -t') + P j^' ds^BBit - s)j-^Gpp{s - t') (117) 

where we assume t > t'. Putting \t'(t — t') back into Eq. flll2p and setting Gbp = (due to 
t > t'), we then have the kinetic equation 

^Blit)Gppit - t') = -f3J:BB{0)Gpp{t -t') + /3 j'^ ds^BBit - s)^Gpp{s - t'). (118) 

This is the same form derived in SDENE, but with the division of the vertices more fully 
defined. 

To continue, we need explicit forms for our local and memory function contributions. 
Collecting terms, we find 



-1 



(5Dpq^ [dt 



V{q) 



/3p2 J [27rY 



V\k)S{k). 



(119) 



The 55-contribution to the memory function is a bit more complex and we address it next. 



B. Memory function 



Let us set up a Fourier transform for the loop pieces which make up the memory function 
in the form, 



MB), 



duj 



^BB v.,t) = - / ^e-*0(^^) 



/3-2ini. 



U3 + — u — ir] 



(120) 
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where AB is short for GG in the case of the collective contribution and GG^^'^ for the single 
particle contribution. 

We first concentrate on the argument 

[1 + iiusK^s + oo^Ki^)? 



Jbb = Im- 



uj^ + — io — irj 



If we change variables such that 



and 



X = UJ3K13 + LJ4Kii 



(121) 



;i22) 



U^U — {us +UJ4), 



we have 



BB 



Im 



1 + 2ix — x'^ 



— — Im 



—U — 17] 

—irj + irjx^ + 2ixu + u + 2xrj — ux"^ 



+ rf 



77(1-^2) 



2xu 



V? -\-rf -\- rf 

Returning to our Fourier transform, we have 



77(1 - x^) 



2xu 



2tx 



,2 I ^2 



-1(^3+0)4)* 



1^1 — x^) + 



(123) 



(124) 



(125) 



At this point, we have now exactly the exponential weight required (cxp(— i(cj3 + u:4)t) 
to factorize the problem. We have (upon substituting x back in) the explicit result 

1 f d'^ks d'^ki 



^BB 



2^2p4y (2ny{2n)d 



{2TiY5{q^ - fcs - k^) 



^^g-i(c.3+a;4)t 



27r 27r 



X 



1 + 2i{ujzKi2, + oJaKia) - {oJzKiz + uJAKnf 



A{k3,u}3)B{kA,u}^) 



1 



d'^ks d'^ki 
2/32p4j (2^(2^ 

X 



{2T:)%q^ -ks- ki) j dt[S{ti - t[) 



l-2K^s^-2K^^— + K^ 
at I at\ 



d_d_ 



(126) 
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X 



where we have introduced an auxihary variable t'^ to help show that the derivatives act only 
on specific terms. 

Moving from generic to our specific loop contributions, we have 

2G,,{h, h)Gfp ih, t\) + Gpp{h, h)G,p{h, t[) (127) 

We see that we generate time derivatives because the three-point vertices are frequency 
dependent. 

In addition to the full time-dependent form, we also need the t = contribution^. In 
the limit of small time, the derivative terms vanish and we find the simple result 

1 f d'^k ~. 



Using our results for the static structure factor to second order given in Eq.f llOSp . we can 
rewrite the second term, giving us 

C. Final form 

Inserting these results into the kinetic equation, we see that a number of terms will cancel 
yielding 

^^G,,{q,t-t') = -Dq'pS-\q)G,,{q,t-t') 

-Dq^ J^^ dsP^pEBB{q,t- s)j-Gpp{q,s -t'). (130) 
Performing a simple shift in time, we can rewrite this as 
§^G,,{q,t) = -Dq'pS-\q)G,,iq,t) 

-Dq^ ds/3^pEBB{q,t-s)^Gpp{q,s) (131) 
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This equation is of the memory function form where we now have a field-theoretic pre- 
scription for the determination of T^BBili^ ~ s). The static part of the memory function 
yields a term proportional to the inverse static structure factor. The dynamic part of the 
memory function is just the BB matrix element of the loop contributions, 7^^''°°^'' + K^if ■ 

(2) 

In SDENE we showed explicitly that Kl- satisfies a FDR and we will show elsewhere that 
(2,«oop) gxplicitly satisfies a FDR as well. 



V. CONCLUSION 



We have shown here that, in the case where one is in thermal equilibrium, the density- 
density correlation function satisfies a kinetic equation of the same form as in MCT-"— . The 
interesting point is that we can explore the corrections of the relevant memory function, T^bb- 
Since T^bb comes from a detailed microscopic derivation, there are several features which 
differ from the conventional mode coupling analysis. At second-order in perturbation theory, 
we have a structure 




4'1(12) = 2 (132) 

where one has a one-loop structure where the three-point vertices have a frequency depen- 
dence. There are two pairs of effective propagators. One contribution is a product of Gij 
propagators where 

~ '2^^^ik ^ kiG ij + GikOktCf^) (133) 
and the other contribution is from the product of Gij and Gf^ where 

Gij = G^^ o'^x 

Gxyay^G%K (134) 

Thus, the microscopic theory is more involved than MOT. In Appendix |Dl we show that 
Gij and Gij themselves satisfy a FDR. 
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At the next order in perturbation theory, one generates two-loop structures such as 



4'kl2) = (135) 

In a companion paper-, we will look at the long-time kinetics generated by the kinetic 
equation using analytic techniques. The main result is that one finds, as in MCT, that the 
late time decay shows two power-law regimes governed by exponents a and h. We show in 
the current case that a and h satisfy the relation 

ni^.A = IIl±^ (136) 

r(i-2a) r(i + 26) ^ ' 

where A is a parameter determined in the model. 

In a second companion paper-, we look at the numerical evaluation of the second-order 
kinetic equation derived here. 
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Appendix A: Summary of the zeroth-order three-point vertex functions 

Let us begin with a few definitions. First, we have the commonly occurring wavenumber 
combinations 

Kij = Dqi ■ qj, (Al) 
= Kii = Dql (A2) 



and 
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We also have the important combination 

G^i = -ioji + Ki. (A4) 
The two-point vertex functions were first derived in SDENE and are 

7<°'(12) = 0, (A5) 



and 



7^(12) = --^5(1 + 2). (A8) 

The three-point vertex functions were first derived in Ref. il3i and are 

7Si(123) = 0, (A9) 
7it(123) = -^^[K,2G,''* + K,sG^''*]5il + 2 + 3) 

= -^[l-«Ei]5(l + 2 + 3), (AlO) 

PPo 

7i°ip(123) = -2-^Ku6il + 2 + 3), (All) 
P Po 



and 



where 



7^^(123) = (A12) 



El = UJ2K12 + CO3K13 (A13) 



and where the other vertices are easily constructed by symmetry. 



13 and are the sum of a reducible 



The four-point vertex functions were also derived in Ref. 
and irreducible contribution, 

(0) _ M (OJ) 

where 

(0,R) _ (0)^(0) (0) (0) ^(0) (0) (0)^(0) (0) . , , .^ 
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The full (amputated) results are 



.(0) = 



-(0) 
ISBBp 



and 



where 



Mt 

Mi2 



(A16) 
(A17) 

(A18) 
(A19) 

(A20) 

Mi2 + Mi3 + Mi4 + M23 + M24 + M34, (A21) 
^PoG^^G-,^K^2[Gl+,{Gl + Gl) + G3+4(G3 + G^)]5{1 + 2 + 3 + 4), (A22) 



IbIp = 2r'G^'*G^^*G-,'*[NT + M12 + Mi3 + M^] 



and 



Nt = |Po{-f^l2-f^34[G'3+4 + Gz+i\ + i^l3-f^24[G2+4 + (^2+4] 

+iri4/r23[G;+3 + Gs+sll^ll + 2 + 3 + 4), 

The irreducible pieces are given by 

725(1234) = 0, 
JO,/) 



7^^(1234) = 0, 



7i°iL(1234) = 0, 



7K,(1234) = 



2 



(5(1 + 2 + 3 + 4) 
5(1 + 2 + 3 + 4) 



-f^l,l+4-f^2,l+4'«l+4 + -f^l,l+3-f^2,l+3'«l+3 
Ql+4'*l+4 + Ql+3'*l+3 



(A23) 

(A24) 
(A25) 
(A26) 

(A27) 



and 



7S?p(1234) = -^5(1 + 2 + 3 + 4) 2 - ia;2(Xi,i+4^2,i+4«i+4 + ^i,i+3^2,i+3«i+3) 

-ia;3(^l,l+2-^3,l+2«l+2 + -^l,l+4-^3,l+4«l+4) 
-iCi;4(^l,H_2-S'4^H-2/«l+2 + -S'l,l+3-S'4,l+3'«l+3) 

|"2-^«,^l (A28) 
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where 



<5l+3 = ^1,1+3^2,1+3, (A29) 

gi+4 = i?l,l+4i^2,l+4, (A30) 

a2 = -ft'l+3,1-^1+3,2 + -^^1+4,1-^1+4,2, (A31) 

"3 = -^1+2,1-^1+2,3 + -S'l+4,l-f^l+4,3 (A32) 



and 



a4 = Ki+2,iKi+2A + -f^i+3,1-^1+3,4- (A33) 

Appendix B: Reduction identities 

A number of "reduction identities" are used in this work. These are part of a larger 



collection which is discussed in Ref. 
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The proofs of these relations are essentially brute 
force demonstrations and it will be sufficient to sketch the basics here. The reduction 
identities fall into two types. 

First, zeroth-order cumulants of a particular number of fields can be reduced to a cumu- 
lant of a smaller number of fields when the argument of the -B-field is zero. In this work, we 
particularly make use of two relations, 

G;g(10) = (Bl) 

and 

G%{1200)=(3'G[f{12), (B2) 

which can be verified by explicitly writing out each cumulant and setting the relevant vari- 
ables to zero. (For convergence, we must set the frequencies to zero first, then the wavenum- 
bers.) The calculation is tedious, but straightforward and so is omitted. 

Second, zeroth-order vertex functions of a particular number of fields can be reduced to 
a vertex function of a smaller number of fields when the argument of the p-field is zero, for 
example with the common term 

= -7Sni2)/po. (B3) 
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Again, this is a straightforward exercise. 

Prom these two simple facts, we can derive other useful identities. When one convolves 
a quantity with a one-point cumulant (either Gi{l) or Gf\l)), the implicit constraining 
delta function will usually cause one or more fields in the product to vanish. For example, 
a combination which appears several times in our work is 

7SUl23)G(°)(3) = 7SUl20)G(°)(0) 
= (-7S(12)/po)po 

= -7Sni2). (B4) 

As another example, consider the first-order contribution to the equation of state (which 
appears again as a component of the second-order contribution) 

= <i(12)a5,(23)G,(3) 
= G(°i(12)aBp(20)G',(0) 
= G(°i(10)aBp(00)Gp(0) 
= G%{W)V{0)p 
= -^6-^^(0)^ 

= -/3popV{0)6{l). (B5) 



Appendix C: Static contribution from Gpp(l) 



Let us look at the integral 

du!i ^ ^ ^ f dui 

27 



<i(l)^(l)GpB(l)V^(l)G(°)(l) 



+GfJ{l)V{l)GBp{l)V{l)G 



(0)/ 



-rG';i{l)V{l)Gpp{l)V{l)G%l{l) 
V\q,) j ^[G(°;(l)(G(°i(l)G,^(l) + G^,(l)Ggi(l)) 

+Gi°i(l)Gpp(l)Ggi(l) 



(CI) 
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Recall the forms of the zeroth-order cumulants 



OJl — IKi 



G<°i(l) = 



and 



where Ki — Dq\. 

Using these, we have 



(C2) 
(C3) 

(C4) 



27r 



doji 



+ IKi 



Gpp(l) 



^1 — 



+/3 



X 



Wi + IKi U\ — IKi 

—IKl 



IKi 



GpB(l) + —Gbp{1) 



(C5) 



For the GpB term, we close the contour integral in the upper half plane while for the Gsp 
term, we close in the lower half plane. This gives 



/ 



27r 



27ri , .^—iK^in^ I ■ \ ~27ri,., %kx . . x 



271 



27r 



-2iKi 



271 \uf + nf 



Gpp(l) 



- ^^GpB{qi,iKi)GpB(qi, -iKi) 



Gpp(l) 



271 \uf + 

Using the FDR to express Gbp and GpB in terms of Gpp, wc have the final result 



(C6) 



/ t^^Pp(l) - PlP'y\<l^) I ^^Gpp{q„., 
+Gpp{qi,u;i 
= plP'V\q,) 



IKi — Ui —IKi — Ui 



1^1 

2 , 9^pp[(ll^^l) 



27r u}( + nf 
= plP'V^qi) J ^Gpp{q,,UJ,). 



The integral over Gpp is simply the static structure factor and we have 



/ 



pt^'V\q^)S{q,). 



(C7) 



(C8) 
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Appendix D: Fluctuation-dissipation relations for G and G 

In this appendix we prove that the dressed propagators G and G individually satisfy the 
same fluctuation-dissipation relation that G satisfies. These results hold at all orders of 
perturbation theory. 

1. G fluctuation-dissipation symmetry 

Recall the form of Gij given by 

^ij = ^{Gfx^xyGyj + Gi^a^yG^y^) (Dl) 

Exphcitly, this yields 

G,B = \{G%VG,B + G,bVG%) = Gf^VG.B. (D2) 
Gb, = \{G%IvGb, + Gb,VG%1) = G%VGb„ (D3) 



and 



G,, = liG%VG,, + G^JvGb, + G.bVG^J + G.^VG^^l)- m 
If we write out our contributing terms as real and imaginary components, 

Gf^^Ro + ilo (D5) 

GpB ^ R + il (D6) 

Gf^ = Ro-tIo (D7) 

Gbp ^ R-il (D8) 

GfJ = ~Io (D9) 

Gpp = --/ (DIO) 

then we have for the imaginary part of Gps, 

ImGpB — —{GpB — Gbp) 

V 
~ 2i 

^V{IRo + IoR). (Dll) 



{Ro + Uo){R + il) - {Ro - ih){R - H) 



Looking next at Gpp, we have 



G 



pp 



2" 



(i^o + ^/o)( --/) +(i?-^/oJ 



00 



+{R + il)(^- ^/o^ + {Ro - tlo) (^-^I 
2V 

^ {2RoI + 2RIo). 

CO 

This imphes, then, the normal fluctuation- dissipation relation, 



Gpp = ImG 

CO 



pB- 



2. G fluctuation-dissipation symmetry 



We may repeat the same procedure for Gij given by 



Explicitly, we have 



and 



pB — '^pB^ '^pBV <^pBj 

CI c^^\rc T//^('-') 

^Bp — U^pV UBpV <^Bp^ 



G,, = G^^JVGnpVG^^l + G'^IvG,,VG'-^1 + G%VG,bVG^^^. 



Using the same decomposition into real and imaginary components, we have 



ImG^B = i ( G 



fpB ~ ^Bp 



2i 

- Yl 

~ 2i 

= V^[I{Rl-I^,) + 2RRoIo] 



{Ro + iIo)iR + ^/)(i^o + ^/o) - (Ro - iIo)iR - ^/)(i^o - 



and 



Gpp = --W^Ro + iIo)(R + il) - V\Rl + - -hV\R - iI){Ro - 



UJ 



{Rl-I^)I + 2IoRoR 



Therefore, we again get the expected fluctuation-dissipation relation, 

2. 



Gpp = ImG 

CO 



pB- 
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Appendix E: Equation of State 

We have identified the equation of state to second order in the pseudo-potential as 



Po = Po(p) = pexp 



(El) 



1. Conventional Form 



To connect this to a more conventional form for the equation of state recall that we are 
in the grand canonical ensemble and 

l^p^ ^ e^i^ (E2) 

where ^ is the chemical potential and I is some microscopic length. We then have the 
thermodynamic identity 

where P is the pressure. Starting with 

/3/. = ln(poO, (E4) 

we have 

d{^P) _ _a(/3//) _pdp, 
dp dp po dp ' 

At first order, everything can be cleanly worked out. Starting with 

Po = pexp[l/(0)] (E6) 

we have 

^ = exp[V^(0)]+pexp[V^(0)]^^ 
= exp[y(0)] + V(0) exp[\/(0)] 



Po 
P 



1 + 1/(0) J exp[1/(0)] 

1 + y (0) ) . (E7) 



Putting this into Eq. (lE5p . we find 



op Po P \ J 

Therefore, 

/3P = p + i/3F(0)pl 

Clearly, we find the ideal gas law and the first order correction. 
If we write more generally that 



we have 



exp[Ty] +pexp[Ty]-^^ 



dp 



dp 



exp[M/] 1^1 + p 
P\ dp 



dp 



which yields. 



dipP) dW 
^ ' = 1+p- 



dp 



dp 



In perturbation theory, we have 



W = V{0) - 



1 r d'^k 



2 J i27iY 



l3'pV\k)S{k) + 



which gives 



dW W I f d'^k 



dp p 2 J (27r) 



P'-pV\k)^^ + ... 
dp 



Recalling that 



S{k) 



1 - pc{k) ' 



we have 



dS{k) ^,„^,d{pc{k)) 



dp 



S\k)- 



dp 
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which finally gives 
dp 

or, integrating, 

fiP = 



S{k) + S\k) 



2.,^^d{pc{k)) 



dp 



d^k 



(5'V\k) 



pS{k) + pS\k)c{k) 



I" 2'^ ' 2 J {27rY 



{k)S\k). 



(E17) 



(E18) 



2. Comparison with Carnahan-Starling Form 



We can go one step further and compare our results with the Carnahan-Starling equation 
of state. This form is an approximate, but quite accurate, equation of state valid for hard 
spheres. 

Recalling the Carnahan-Starling forn>i^ii^, we have 

l3P _ 1 + 7] + 7]^ -7]^ 

p (i-v)^ 

Let us begin by taking the derivative of this with respect to p, 

d{PP) d f _1 + 7] + 7]^ — 7]'^ 

dp 



(E19) 



d-p^P- 



{1-7]Y 



1 + 7] + 7]"^ — 7]^ d f 1 + 7] + 7]^ — r]^ \ dr] 



(1 — r]Y dr]\ (1 — 7]Y J dp 

1 + ?7 + - f {1 - 7]){1 + 27] - 37]'^) + 3(1 + 7] + 7]'^ - 7]^^ 



1 + 4?7 + 47/2 - 4:7]^ + 7]"" 



{i-7]y 



(E20) 



{1-7]Y 

This result is the left-hand side of the thermodynamic relation in Eq. flESp . Let us now 
rewrite the right-hand side. 
We have 

p dpo p 



Po dp po 



'-^wM 9Wiv) d7] ^ ^^(^; 
d7] dp 



dW{7]) ^ 
^ d7] 



(E21) 



where we again use the general form for po given by Eq. flElOl) . 
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Setting the left and right halves equal, we have 

1 + 4r7 + 477^ - 4?7^ + ri^ dW(r]) 

— = v- 



+ 1. 



(1 — T])^ ' dr] 

which we may rearrange as 

dW{ri) 1 + ir] + Ar]"^ - 4r]^ + 7]^ - {1 - T])^ 2(4 - r/) 



drj 

Integrating, we find 



W{r]) = / dx 



2(4 -x) _8r]- 'drf + Sr^^ 
' [l-xY ~ (1-77)3 



or, returning to the full form for po. 



Po = pexp 



8?7 - 'drf + ?,7f 



(E22) 



(E23) 



(E24) 



(E25) 



(1-7^)3 

We now have an independent measure for the quality of our equation of state results. 
As one self-consistently solves for the pseudo-potential, we may compare our perturbative 
result Eq.f lEip to this result. 
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